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Abstract 

The adsorption of a collagen fragment on both a hydrophobic, hydrogen-terminated and a hy- 
drophilic, natively oxidised Si surface is investigated using all-atom molecular dynamics. While 
favourable direct protein-surface interactions via localised contact points characterise adhesion to 
the hydrophilic surface, evenly spread surface/molecule contacts and stabilisation of the helical 
structure occurs upon adsorption on the hydrophobic surface. In the latter case, we find that 
adhesion is accompanied by a mutual fit between the hydrophilic/hydrophobic pattern within the 
protein and the layered water structure at the solid/liquid interface, which may provide an addi- 
tional driving force to the classic hydrophobic effect. 
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The issue of protein adsorption onto solid surfaces in a wet environment is of great 
importance for the design and functionalisation of materials in the contexts of biomedical 
implants, biosensors, antifouling surfaces, pharmaceutical packaging, and many others 
jj] . One important case is the adsorption of extracellular matrix proteins such as collagen on 
the surfaces of implanted Si-based microelectromechanical systems j^J]. Control over protein 
adhesion may be achieved by manipulation of surface properties, such as isoelectric point, 



functional group termination, topography and hydrophobicity 
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However, the atomistic details of the interactions mediated by these effects are still unclear 
and are the object of intense experimental and theoretical investigation yy y y. 



In a previous work [17(, we have shown that liquid water close to a natively oxidised Si 
surface presents alternating layers of larger and smaller density with respect to the bulk. 
These structural features directly influence the strength of mutual adhesion between pairs 
of Si wafers at room temperature under humid conditions 17|. In a similar way, the ad- 
sorption and adhesion of biomolecules to solid surfaces may be significantly determined by 



the structure of water layers in the surface proximity 18 
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21]. In particular, the dis- 



ruption of the adsorbed water layer by individual amino acids has been shown to have both 



enthalpic and entropic effects on the free energy of the bound system [22j. In this work we 
show through all-atom molecular dynamics simulations (MD) 23| that a complex oligopep- 
tide interacts with a Si surface via a mutual fit between its intrinsic hydrophilic/hydrophobic 
pattern and the structured water layers close to the surface. This may provide a driving force 
for ad sorp tion at the solid/liquid interface which acts in addition to the classic hydrophobic 



effect 
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24]. 



Specifically, we use MD to investigate the adsorption behaviour of the NC1(84-116) do- 
main of collagen XIV (PDBID: lb9q), on both a hydrophilic (i.e. natively oxidised 25j) 
and a hydrophobic (i.e. hydrogen-terminated) Si(100) surface. The a-helical protein confor- 
mation used as the starting point for our simulations has been found by circular dichroism 
measurements to be stabilised by interaction with hydrophobic surfaces 26|. The folded 
collagen structure presents one predominantly hydrophobic side and a series of hydrophilic 
basic residues (lysine) clustered on the opposite side. These basic residues form a specific 
binding site through which collagen XIV interacts with heparin molecules 26|. Here, we 
consider adsorption of the a-helix on solid surfaces with the heparin binding lysine residues 
pointing towards the solvent, away from the surface. 
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In our simulations, water molecules are treated using the TIP3P force field 27| and 
protein-water and protein-protein interactions are described by the Amber parm99 biomolec- 
ular force field j^. The water solvent and the protein interact with the surface via Coulomb 
and Lennard- Jones (LJ) interactions. The atomic charges and interatomic interactions for 
the natively oxidised Si surface are taken from Ref. [nj]. The LJ energies and radii oi of 
the Si and O surface atoms are accurately tuned to reproduce the correct heat of immersion 

n 

of hydrated silica surfaces [17jj. I n t ne case °f the hydrophobic H-terminated Si(100) sur- 
face, terminal H atoms are assigned small ESP charges of -0.08 e (as computed by density 
functional theory), which are compensated by charges of +0.08 e on their neighbouring Si 
atoms. For the H atoms of the hydrophobic surface, we use the same LJ parameters as 
for hydrophobic H atoms of hydrocarbon chains. The surface-protein LJ interactions are 
obtained by combining the surface and biomolecular force field parameters using the stan- 
dard Lorentz-Berthelot combination rules 29| . We have performed both shorter simulations 



allowing all atoms of the surfaces to move and longer simulations where the surfaces remain 
frozen throughout, except for hydroxyl H atoms on the natively oxidised surface, which are 
free to rotate as described in Ref. 17| . The water structure at the rigid natively oxidised 
surface is nearly identical to that obtained for the unconstrained surface, as also reported 
in [17|. 



Two simulations are started initially placing the collagen helix a few A above two ran- 
domly chosen sites of the natively oxidised Si surface 23J . Close contact between the surface 
and collagen is achieved by performing a 2000 step conjugate gradients minimisation in vac- 
uum. The entire simulation cells are then filled with water molecules and the cell heights are 
adjusted to obtain the equilibrium TIP3P water density far from the surfaces. The solvated 
systems are slowly heated to 300 K applying weak harmonic constraints to the protein and, 
finally, constant-temperature production MD runs are performed for 24 ns. 

The collagen fragment is found to be highly mobile on the hydrophilic surface 
(Fig. CD (top)). In the first simulation, the molecular centre of mass initially moves away 
from the surface, while localised contact (within 3.4 A of the surface) is maintained via 
residues 108-112, resulting in an end-on bound conformation. Later a second point of con- 
tact via residues 91-93 is formed. The position of the centre of mass remains stable at 
approximately 9 A over the final 10 ns of the simulation, and the root mean square displace- 
ment of the backbone atoms of the protein converges to within 1 A. In the second simulation, 
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two points of contact are maintained throughout most of the 24 ns simulation, again via the 
molecular termini (roughly localised in the regions of residues 84-88 and 104-112). Inter- 
estingly, residues Leul07 and Arglll remain in close contact with the surface throughout 
both simulations via the hydrophobic CH n part of their side chain groups. The total energy 
of adsorption averaged over the final 10 ns of the two simulations is —22 ± 7 kcal/mol with 
respect to a reference unbound simulation in which collagen is separated from the surface 
slab by at least 10 A. 

A different behaviour is observed in two similar simulations (each lasting 14 ns) performed 
starting with collagen in the same conformation over two different starting positions close 
to the hydrophobic, H-terminated Si surface. The fragment remains adsorbed to the surface 
and maintains its helical structure throughout both simulations (Fig. [T] (centre)), again with 
an average energy of adsorption of — 22±3 kcal/mol with respect to a corresponding unbound 
system. Contrary to the case of the hydrophilic surface, the residues in direct contact with 
the surface are evenly spread through the central part of the protein (residues 91-108), while 
the molecular termini (84-90 and 109-116) are less tightly bound. 

As visible in the Ramachandran plot computed over the final 10 ns of each simulation 
(Fig. [2]), the backbone dihedral angles of collagen on the hydrophobic surface are tightly 
clustered around the ideal a-helix values (-65°, -35°). A time-averaged calculation of the 



301 ] reveals that 33 % of the 



a-helix content of the protein using the DSSP algorithm 
residues have an a- helical structure. This is mainly concentrated in the central protein 
region of residues 93-112, in qualitative agreement with a structural assignment based on 
NMR spectroscopy (^fj. For comparison, in a separate 14 ns simulation of collagen in bulk 
water, there is a much wider spread of dihedral angles, with an increase of their occurrence 
in the /3-sheet region. With respect to the adsorbed case, the average helical content falls 
to 18 %, and the average number of intramolecular hydrogen bonds decreases from 19 to 
17 (as computed with a distance cut-off of 3.4 A and an angular cut-off of 140°). On the 
hydrophilic surface, the degree of unfolding is intermediate between the hydrophobic and 
unbound structures (Fig. [2]) and 20 internal hydrogen bonds are maintained on average. 

Although the degree of a-helical conformation in bulk water may be over-estimated due 
to our helical starting structure, the uncoiling observed in the time scale of these simula- 



tions appears to be significant. In Ref. 



261 ]. it is contended that unfavourable electrostatic 



interactions between arginyl and lysyl side chains in the region of the heparin binding site 
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contribute to helix destabilisation which causes the collagen domain to unfold in water. 
Indeed, by computing the pairwise electrostatic interactions involving only the primary 
heparin binding region (residues 97-105), we find that the net electrostatic energy is approx- 
imately 30 kcal/mol less favourable at the hydrophobic surface than in bulk water, where 
the Coulomb repulsion between charged residues is lowered as a consequence of the helical 
unfolding. 

These results are consistent with the experimental observations that the a-helical form 
of the NCI collagen XIV fragment adsorbs stably on hydrophobic organic surfaces 26]. 
Here we observe that collagen remains adsorbed at the hydrophobic Si-H surface for a total 
of 28 ns with a favourable net binding energy. On the one hand, displacements of water 
molecules close to non-polar amino acids following direct interactions with the surface, to- 
gether with the associated re- arrangement of the solvent hydrogen-bond network, are indeed 
known to drive adsorption and represent a fundamental part of the so-called hydrophobic at- 



traction 
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24j . On the other hand, more subtle structural effects arising from the presence 



of an abrupt solid/liquid interface combined with a favourable arrangement of the adsorbed 
protein in an ordered hydrophobic/hydrophilic pattern may also contribute to stabilize the 
adsorption, as we shall see below. 

The effects of the molecular water structure on the biomolecular adsorption are studied by 
analysing the density profiles of water and protein H atoms along the direction perpendicular 
to the surface. In the case of water, this is computed by counting the average number of H 
atoms present in thin planar slices of width 0.1 A in separate 1 ns simulations in the absence 
of collagen. The disruption of the bulk water hydrogen bond network by the two surfaces 
results in a re-ordering of interfacial water molecules to maximise their interactions with 
neighbouring molecules and the surface. This is manifested by the structuring of water into 
layers of higher and lower density with respect to the bulk, as found previously in numerous 



MD simulations 
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31 



32), 



In the case of the hydrophobic surface (Fig. [3]^a), black line), the first peak in the water 
density is separated from the surface H atoms (grey lines) by 2.7 A, consistent with weak 
surface-water interactions of the Lennard- Jones type [33}]. Two other maxima are visible at 
5.9 and 9.1 A far from the surface. This surface-induced water layering is preserved in the 
presence of adsorbed collagen, as evident from the density profile of the 150 water molecules 
belonging to the first two solvation shells of the protein fragment (dotted line). The first 
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three peaks are centred at 2.8, 6.5, and 9.7 A, only slightly shifted from those observed 
in the absence of the protein. We now compute in the same way the density of H atoms 
belonging to collagen, averaged during the MD trajectories described above, discriminating 
hydrophobic H atoms (bonded to C atoms) from hydrophilic H atoms (bonded to O, N or 
S atoms). The H atoms of the protein are clearly structured in alternating regions of higher 
and lower density, revealing an intrinsic ordering of hydrophobic/hydrophilic residues in the 
adsorbed helical structure (Fig. E^a)). Interestingly, the hydrophobic/hydrophilic pattern 
within the adsorbed protein is commensurate with the water layering oscillations in the 
surface proximity. Namely, the first water peak is located between the first and second 
small peaks of the hydrophobic H atoms, centred at 2.2 and 3.9 A. The latter peak coincides 
approximately with the first trough in the water density (located at 4.4 A). Correspondingly, 
the large trough in the hydrophobic H density at 6.2 A coincides with the second water 
density peak. Further, the first peak in the hydrophilic H density is located at 2.8 A, in 
perfect correspondence with the large first water peak. The broad main peak overlaps with 
the second water peak and the third peak (corresponding to the cluster of lysine residues) 
points away from the surface/protein interface into the bulk solution 34]. 

In the case of the natively oxidised surface, water molecules penetrate into the oxide 
layer and the main density peak is centred on the surface hydroxyl groups (Fig. [3](b)). 
In contrast to the clear structuring observed at the hydrophobic surface, both hydrophilic 
and hydrophobic H atoms of the protein show relatively spread and less structured density 
profiles, and the matching between the hydrophilic/hydrophobic patterns in the molecule 
and the water layer structure is barely visible. Surprisingly, while hydrophilic H atoms of 
the protein are shielded from the surface by the first high density water peak, the H atoms 
of hydrophobic groups are the closest to the surface. This is consistent with our previous 
finding that predominantly hydrophobic side chain groups remain in close surface contact 
(see above), bound via van der Waals forces to bridging oxygen atoms of the oxide network, 
which are hydrophobic compared to the surface silanol groups 

Despite the evident ordering of water molecules at the surface/protein interfaces, their 
mobility within the protein hydration shells remains of the same order of magnitude before 
and after adsorption. This is shown by the mean residence times of the first shell water 
molecules surrounding the wet surfaces, the adsorbed proteins and the protein in bulk water, 
calculated with the method introduced in Ref. 35] (Table H]). As expected, the residence 
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times of water molecules in the surface proximity, especially for the hydrophilic case, are 
higher than the value of 7 ps obtained for bulk water. However, the residence times of water 
molecules around collagen (80 ps for the dissolved protein, as found for similar systems 351 ]) 
are only slightly increased when the protein is adsorbed at the two surfaces. 

To summarise, collagen remains bound to the oxidised Si surface via a few terminal 
residues, the most stable direct interactions with the surface being between hydrophobic 
side chain groups and hydrophobic surface sites. Individual hydrophilic residues penetrate 
into the first water layer to form direct hydrogen bonds with the surface for only 20-30 % 
of the total simulation time. This is consistent with a suggested adsorption model in which 
proteins weakly bound to hydrophilic surfaces leave the surface-bound hydration layer almost 
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In contrast, 



12|. 



intact and do not affect the surface tension at the solid/liquid interface 
proteins tend to adsorb with larger contact surface areas to hydrophobic surfaces [5|, [111, 
In particular, the NCI domain of collagen XIV has been observed to adsorb stably in an 
a-helical structure on a hydrophobic organic layer [26j, as supported by our simulations on 
the H-terminated Si surface. 

We note that an accurate analysis of the individual enthalpy and entropy contributions 
leading to the observed adsorption behaviour is made difficult by the large oscillation of each 
separate energy term and the intermittent exchange of counter-ions between surface, protein 
and bulk water. Extensions of this work will be necessary to calculate the full potential of 
mean force along possible adsorption/desorption paths at different temperatures In 
our analysis so far, we have focused on the structural details governing collagen adsorp- 
tion on silicon and provided evidence for a strict interplay between the solvent structure 
at the solid/liquid interface and the intrinsic hydrophobic/hydrophilic pattern within the 
adsorbate. Whether this effect is specific to systems that are intrinsically chemically com- 
mensurate with the water layering, or whether flexible biopolymers may indeed adjust their 

in a process rem- 



371] ) is a question 



internal structure to fit the solvent ordering at the solid/liquid interface 
iniscent of the "induced fit" mechanism for supramolecular organisation 
that will require further combined theoretical and experimental efforts. 
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Figure Captions 

FIG. 1: MD trajectories of collagen on a natively oxidised (top) and a H-terminated (centre) 
Si(100) surface and in a separate simulation in bulk water, in the absence of any surface (bottom). 
Initial (left), intermediate (centre) and final (right) snapshots are pictured. Collagen consists of 
34 amino acid residues and its predominantly hydrophobic surface is initially oriented towards the 
surfaces. Water molecules are omitted for clarity. 

FIG. 2: Ramachandran plots showing the distributions of the dihedral angles, <j> and ift (degrees), 
for collagen at the hydrophobic surface (left), hydrophilic surface (centre) and in bulk water (right). 

FIG. 3: Density profiles of hydrophobic (green) and hydrophilic (blue) H atoms of collagen at (a) 
a H-terminated and (b) a natively oxidised Si surface. Density profiles of water H atoms (black, 
solid line) are obtained from separate 1 ns simulations. Also shown are the H atoms of water 
molecules belonging to the first two solvation shells of collagen (black, dotted line) and H atoms 
of the surfaces (grey). 

FIG. 4: Water mean residence times (ps) about the protein adsorbed at the hydrophobic and 
hydrophilic surfaces and about the protein in bulk water. Also shown are the residence times of 
water molecules at the two surfaces in the absence of collagen. The two values quoted at hydrophilic 
surfaces are for hydroxyl groups and bridging oxygen atoms respectively. Snapshots are analysed 
every 2 ps in sliding windows of 400 ps in the final 10 ns of each simulation. The mean residence 
times are obtained by fitting a stretched exponential to the autocorrelation functions of the boolean 
variable that indicates whether or not a given water molecule is within a shell of width 3.5 A around 
the protein or surface. 
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Figure [U D. J. Cole et al. 
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Figure EJ D. J. Cole et al. 
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Figure El D. J. Cole et al. 
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Figure H D. J. Cole et al. 
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